Terahertz large-area unidirectional surface magnetoplasmon and its applications

Unidirectional surface waves in nonreciprocal plasmonic platforms with nonlocal effects have been a topic of significant interest and some controversy. In this study, we present a scheme to achieve unidirectional surface magnetoplasmons (USMPs) with large modal areas at terahertz frequencies. Such large-area USMPs (LUSMPs) exist in a metal-UENZ (uniaxial-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon $$\end{document}ε-near-zero)-Si-InSb structure under external magnetic field, where the effect of nonlocality is included. The field of the LUSMP extends almost uniformly in the UENZ layer with a thickness of wavelength scale, thus its modal size can be represented by the UENZ-layer thickness. Due to the modal energy primarily distributed in the thick UENZ layer, the nonlocality-induced leakage of the LUSMP is significantly reduced by an order of magnitude, compared to previous USMP existing at interface between InSb and opaque isotropic medium. Due to their large modal sizes, such LUSMPs can be efficiently excited by terahertz radiations directly from free space. In addition, LUSMPs offer high degree of freedom for manipulating terahertz waves, such as energy squeezing and trapping. Based on LUSMPs, a terahertz free-space isolator is also developed. Our findings have important implications to the development of innovative plasmonic devices in terahertz regime.


Terahertz large-area unidirectional surface magnetoplasmon and its applications
Qian Shen 1,2 , Jinhua Yan 1 , Yun You 1 , Senpeng Li 1 & Linfang Shen 1* Unidirectional surface waves in nonreciprocal plasmonic platforms with nonlocal effects have been a topic of significant interest and some controversy.In this study, we present a scheme to achieve unidirectional surface magnetoplasmons (USMPs) with large modal areas at terahertz frequencies.Such large-area USMPs (LUSMPs) exist in a metal-UENZ (uniaxial-ε-near-zero)-Si-InSb structure under external magnetic field, where the effect of nonlocality is included.The field of the LUSMP extends almost uniformly in the UENZ layer with a thickness of wavelength scale, thus its modal size can be represented by the UENZ-layer thickness.Due to the modal energy primarily distributed in the thick UENZ layer, the nonlocality-induced leakage of the LUSMP is significantly reduced by an order of magnitude, compared to previous USMP existing at interface between InSb and opaque isotropic medium.Due to their large modal sizes, such LUSMPs can be efficiently excited by terahertz radiations directly from free space.In addition, LUSMPs offer high degree of freedom for manipulating terahertz waves, such as energy squeezing and trapping.Based on LUSMPs, a terahertz free-space isolator is also developed.Our findings have important implications to the development of innovative plasmonic devices in terahertz regime.
Unidirectional propagation of surface magnetoplasmons (SMPs), as a striking phenomenon in nonreciprocal and topological electromagnetics, have attracted growing interest and also caused some controversy.SMPs are nonreciprocal surface plasmon polaritons in plasmonic materials magnetized by external magnetic field (which breaks time-reversal symmetry), and their studies began five decades ago [1][2][3][4] .In recent years, SMPs have received a resurgence of interest due to possibility of realizing truly unidirectional transport channels [5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21] .For nonreciprocal plasmonic platforms, there may exist unidirectional frequency windows, in which SMPs are allowed to propagate only in one direction.Such unidirectional SMPs (USMPs) can even be immune to backscattering at imperfections and defects.Remarkably, when one-way channel is terminated, USMPs might be completely trapped, resulting in an accumulation of electromagnetic energy with high field localization [16][17][18] .Thus, an apparent violation of the time-bandwidth limit is achieved 19 .However, subsequent studies argued that USMPs reported in Ref. 19 are not strictly unidirectional when nonlocal material effects are considered 20 .Recently, true USMPs that are robust to nonlocal effects have been reported 21 .Such USMPs exist at the interface between magnetized plasma and opaque medium and are topologically protected.However, due to the nonlocality, these USMPs become leaky and suffer from serious leakage losses 21 .For a typical example in Ref. 21 , USMP only has a propagation length shorter than two (vacuum) wavelengths even in the absence of material loss, and this would hinder its practical applications.Hence, it is necessary to develop alternative types of USMPs that can completely eliminate or significantly mitigate the impact of nonlocal effects.
As surface waves, the fields of USMPs are generally confined in a narrow boundary.Recently, unidirectional modes with large modal areas have been reported.In the microwave regime, researchers have demonstrated largearea unidirectional modes in pseudospin-field-dependent waveguides and photonic crystals containing Dirac points 22,23 .Large-area USMPs (LUSMPs) have also been achieved in the microwave regime by using a uniaxial µ-near-zero material 24 .For such LUSMPs, most of the modal energy is distributed within the thick µ-near-zero material layer.Therefore, it is naturally desired that LUSMP mode can be extended to terahertz (THz) regime.In recent years, rapid development in terahertz science and technology has been made because of huge potential in imaging, spectroscopy, biomedical sciences, and communications.For those applications, terahertz USMPs have been proposed 7,21 , which are of particular importance for functional devices such as isolators, switches and splitters.Such USMPs can be achieved using semiconductors under external magnetic field, and they may be classified into two types.The first-type USMPs are the conventional SMPs, and exist in the lower bandgap of the magnetized semiconductor 7 , but they have been recently shown to lose their strict unidirectional nature if nonlocal effects are included 20 .The second-type USMPs are true unidirectional modes, existing in the upper bandgap of the magnetized semiconductor, and they are robust to nonlocality because the upper bandgap is topologically nontrivial 21 .However, these robust USMPs still suffer from serious leakage caused by the nonlocality.If LUSMPs could be achieved at terahertz frequencies, the impact of nonlocal effects in plasmonic waveguides would be largely weakened and the modal leakage can be reduced to a very low level.In addition, LUSMPs can also provide high degree of freedom to manipulate terahertz waves.
In this paper, we propose a scheme to achieve terahertz LUSMPs using uniaxial ε-near-zero (UENZ) material 25 .Different from all previous terahertz USMPs, which are tightly localized on the semiconductor surface, the proposed terahertz LUSMPs can possess modal spots with sizes larger than a wavelength.Such terahertz LUSMPs exist in the upper bandgap of the magnetized semiconductor, so they are topologically protected and robust to nonlocality.More importantly, compared to previous (true) terahertz USMPs, the nonlocality-induced leakages of LUSMPs at terahertz frequencies can be significantly reduced by an order of magnitude, because their modal energies are mainly distributed in the thick UENZ layer.Owing to their large modal spots, these LUSMPs can be efficiently excited by terahertz radiation incident from free space.Moreover, based on LUSMPs, we will numerically demonstrate several interesting examples of flexibly manipulating terahertz radiations.

Physical model and dispersion diagram
The proposed waveguide for supporting terahertz LUSMP is a metal-UENZ-Si-semiconductor stratified structure under an external dc magnetic field B 0 , as illustrated in Fig. 1a.We start with analysis of the dispersion for the waveguide in local material model.In this waveguide, the metal is assumed to be perfect electric conductor, which is a valid approximation in terahertz regime.The Si layer with thickness d 1 has the relative permittivity ε d = 11.68 20 , and the UENZ material layer with thickness d 2 has the relative permittivity in the form where ε x is larger than 1, and ε y = ε z is close to zero.The semiconductor with gyroelectric anisotropy induced by the external magnetic field is characterized by a relative permittivity tensor with where ω is angular frequency, ω p is the plasma frequency, ν is the electron scattering frequency, ω c = eB 0 /m * (where e and m * are, respectively, the charge and effective mass of the electron) being the electron cyclotron frequency, and ε ∞ is the high-frequency permittivity.The magnetized semiconductor exhibits two bandgaps for bulk modes with transverse-magnetic (TM) polarization: a lower bandgap below the plasma frequency, which exists even in the absence of external magnetic field, and an upper bandgap above the plasma frequency, which is opened by the external magnetic field.The upper bandgap has nontrivial topological property characterized by a nonzero gap Chern number, and as a result, USMP within it is topologically protected 21 and thus could be robust to nonlocal effects.
The waveguide supports SMPs with TM-polarized field (i.e., H = ŷH y ).The dispersion relation of SMPs in the local material model can be derived analytically from Maxwell's equations combined with proper boundary conditions, and we obtain where k is the propagation constant, are, respectively, the transverse attenuation coefficients of the field in the UENZ material, silicon, and semiconductor.From Eq. ( 3), two asymptotic frequencies can be obtained, at which k → ±∞ , and they are given by ( 1) In this paper, the semiconductor is assumed to be n-doped InSb with ε ∞ = 15.68 and ω p = 4π × 10 12 rad/ s 20,21 , and we set ω c = 0.25ω p , which corresponds an external magnetic field of 0.25 T. Figure 1b shows the dispersion diagram for the waveguide with ε z = 0.001 , ε x = 50 , d 1 = 0.0042 p and d 2 = 2 p (where p = 150 µ m, is the vacuum wavelength for the plasma frequency).There exist two pairs of SMP branches (red lines).The lower pair of SMP branches are located in the lower bandgap of the magnetized InSb, and they have different asymptotic frequencies as well as different upper cutoffs.In the frequency range from the upper cutoff of the right branch to the asymptotic frequency of the left branch, SMP is only allowed to propagate in the backward direction, and this is just a unidirectional frequency window in the lower bandgap.The upper pair of SMP branches have different lower cutoffs, thus creating another unidirectional frequency window in the upper bandgap.For our waveguide, USMP in the upper frequency window exists with small wavenumbers, and it is only allowed to propagate forward.The unidirectional frequency windows of two types are indicated by the yellow shaded areas in Fig. 1b.Here, it should be noted that if the Si layer is removed from our waveguide, the SMP branches in the lower bandgap will extend into the upper bandgap, leading to the disappearance of the upper unidirectional window.Besides, recent studies show that the asymptotic frequencies of SMPs vanish when the nonlocal effect of semiconductor is included, and consequently, the waveguide loses its lower unidirectional window 20,21 .This nonlocal response is a result of the free electron movement during an optical cycle, caused by convection and diffusion.The nonlocal effect can significantly modify the dispersion feature of SMPs, implying that the local model is inadequate.Therefore, in what follows, we focus on USMPs in the upper bandgap of the magnetized InSb.For our waveguide, it is also feasible that the unidirectional window covers the whole upper bandgap, as illustrated in Fig. 1d, where ε x = 300 and d 1 = 0.0012 p , and the other parameters are the same as in Fig. 1b.In the nonlocal model, we adopt a hydrodynamic model of free-electron gas (without considering diffusion effects) to describe the semiconductor.From a microscopic hydrodynamic equation of motion for the free electrons that includes a pressure term, the free-electron current J induced by the field in the semiconductor satisfies the following equation where β is the nonlocal parameter.Here, we take β = 1.07 × 10 6 m/s for InSb at room temperature 20 .The wave propagation in the magnetized semiconductor is governed by Maxwell's equations in the form By simultaneously solving the Maxwell's equations combined with Eq. ( 5), the dispersion equation for TMpolarized bulk modes in the semiconductor can be derived as where k b is the bulk-mode wavenumber.Compared to the local dispersion relation of bulk modes, k b = ± √ ε v k 0 , the nonlocal dispersion Eq. ( 7) implies that there exist two types of bulk modes in the nonlocal semiconductor.The green shaded areas in Fig. 1b depict the bulk-mode dispersion of InSb, which differs from that of the local model.In the nonlocal scenario, the lower bulk mode has no upper cutoff and covers the upper bandgap for the local model, while the upper bulk-mode zone remains similar to that for the local model and still exhibits a high-frequency bandgap.
Therefore, for the waveguide in the nonlocal model, the field of a SMP mode has two bulk-mode components in the semiconductor.In order to solve for the SMP modes in the waveguide with nonlocality, at the interface between the semiconductor and Si, an addition boundary condition of J x = 0 must be imposed, except for the ordinary boundary conditions, which require the continuity of the field components H y and E z .As in the local case, the dispersion relation for the SMP modes in the nonlocal case is analytically derived, and it can be written as with where γ 1 = k 2 − k 2 b1 and γ 2 = k 2 − k 2 b2 ( k b1 and k b2 are the wavenumbers of the two bulk modes in the magnetized semiconductor).The dispersion Eq. ( 8) is numerically solved for the case in Fig. 1b, and the results are also plotted as dashed lines in this figure.In the nonlocal case, the asymptotic frequencies of SMPs in the lower bandgap of the magnetized InSb vanish, resulting in the disappearance of the USMP phenomenon.In the upper bandgap of the magnetized InSb, the dispersion curves of SMPs are identical for both the local and nonlocal cases, and this is more clearly seen in Fig. 1c.Besides, no backward mode emerges within the validity range ( |k| < 150k p ) of the hydrodynamic model.So USMP and its unidirectional window in the upper bandgap are preserved in the presence of nonlocality.This is also the situation for the case in Fig. 1d.However, as the lower bulk mode becomes propagating in the (local) upper bandgap, this USMP become a leaky mode, because its field has a radiation component in the InSb.Consequently, the propagation constant k of USMP becomes a complex number in the nonlocal model.Therefore, it is necessary to further investigate the impact of the nonlocal effect on the modal properties of the USMP.

Mode analysis and transmission simulation
In our waveguide, the USMP has a H-field distribution in the UENZ layer, characterized by α e in Eqs. ( 3) and (8).As ε z close to zero, α e approaches zero as well, allowing for a large modal area when d 2 is large.To demon- strate this, we calculated the H-field profiles of USMP at f = 2.12 THz for both the local ( β = 0 ) and nonlocal β = 1.07 × 10 6 m/s) cases, and the results are plotted in Fig. 2a, where ε z = 0.001 , ε x = 50 , d 1 = 0.63 µ m, and d 2 = 300 µ m.In the UENZ layer, the H-field extends uniformly in the both cases, resulting in a transverse modal size of USMP that can be represented by d 2 .This modal size is nearly two (vacuum) wavelengths, implying that this USMP is a LUSMP mode.By using a UENZ material with ε z = 0.001 and ε x = 50 , the maximal modal size of USMP can reach up to 14.1 p .It should be noted that for the nonlocal model, the H-field amplitude does not decay to zero in the lossless InSb layer.As the propagation constant of LUSMP is smaller than the wavenumber of the additional bulk mode, the LUSMP mode has a radiation component in the semiconductor and thus becomes a leaky mode with a complex propagation constant k = k r + ik i , where k r is the phase constant and k i is the attenuation constant.It should be emphasized that due to the radiation loss, k i is nonzero even if the materials in the system are lossless.Our numerical calculations show that k i is closely dependent on d 2 .Figure 2b depicts the dependence of k i on d 2 for the nonlocal ( β = 1.07 × 10 6 m/s) and lossless ( ν = 0 ) case, where the frequency is f = 2.12 THz, and the other parameters are the same as in Fig. 2a.k i decreases as d 2 increases.As (5) www.nature.com/scientificreports/d 2 increases from 0.1 to 2 p , k i decreases from 0.04 to 0.0035k p .It is worth noting that at the same frequency, k i for the regular USMPs (RUSMPs) 21 is approximately 0.04k p .Obviously, compared to RUSMPs, the radiation loss due to the nonlocality is largely reduced by an order of magnitude for the present LUSMP.For the LUSMP, the modal energy is primarily concentrated in the UENZ layer, effectively weakening the radiation effect in the semiconductor.This makes LUSMPs have a significant advantage over RUSMPs at terahertz frequencies.
In the lossy case ( ν = 0 ), k i has two components: radiation loss due to nonlocal effect and absorption loss due to lossy materials.Figure 2b shows that the total k i of USMP in the nonlocal and lossy cases decreases as d 2 grows.Various values of the loss parameter ν are analyzed, including ν = 0 .Evidently, the absorption loss of USMP can also be effectively reduced by increasing d 2 , because of the same reason for decreasing the radiation loss.To further clarify this, we also present the dependence of k i on d 2 for the local cases with different ν values in Fig. 2c.In the local case, k i is only caused by the absorption loss, and it is indeed found to be decreasing with increasing d 2 .Therefore, our guiding system can effectively reduce both radiation loss and absorption loss by increasing d 2 .We further compare the propagation lengths [ L = 1/(2k i ) ] of the LUSMP and the corresponding RUSMP in Ref.21, which is the only example of true USMP previously reported.Figure 2d shows the propagation lengths of the LUSMP and RUSMP as functions of frequency.For our waveguide, the parameters are as follows: d 1 = 0.63 µ m, d 2 = 300 µ m; ε z = 0.001 , ε x = 50 .At f = 2.12 THz, we find that for ν = 0 , the propagation lengths are L RUSMP = 231 µ m and L LUSMP = 3402 µ m.When ν = 0.002ω p , the propagation lengths become L RUSMP = 166 µ m and L LUSMP = 2883 µ m, implying that the radiation loss is a dominating factor of the LUSMP and RUSMP losses for low-level material loss.However, when ν = 0.02ω p , we find that L RUSMP = 48 µ m and L LUSMP = 1215 µ m, which are close to their values of 77 µ m and 1650 µ m obtained from the local models, implying that the absorption loss dominates in this case.
To confirm the validity of LUSMP in the nonlocal model, we simulated wave transmission in the proposed waveguide ( β = 1.07 × 10 6 m/s, d 2 = 2 p ) using the finite element method (FEM).The simulation employed a magnetic current line source as an excitation, specifically operating at 2.12 THz (i.e., marked by A in Fig. 1c), and the source was placed at the middle of the UENZ layer.A material loss of ν = 0.002ω p was considered.The simulated H-field amplitudes are depicted in Fig. 3a.Notably, the excited wave propagates only in the forward direction, and the H-field is uniformly distributed in the UENZ layer away from the source.To examine the robustness of LUSMP in the nonlocal model, we inserted a defect into the waveguide, and re-simulated the guiding system.This defect (defect I), located in the middle of the UENZ layer, is an air column with a square cross-section of side length 0.5 p .But no backward wave is generated by the defect, and the incident LUSMP completely bypasses it and continues to propagate forward, as illustrated in Fig. 3b.The defect only modifies the local field around it, and behind it the field quickly recovers.This is also clearly displayed in Fig. 3g, where the distributions of magnetic field amplitude along the middle line of the UENZ layer are plotted for both cases with and without defect.To investigate the influence of the defect size and location, we further move the defect down by 0.5 p (defect II), or double its size (defect III).The simulated results are displayed in Fig. 3c,e, respectively.In both the cases, the incident LUSMP bypasses the defect, but the field behind the defect is somewhat weakened, www.nature.com/scientificreports/which is more shown in Fig. 3g.Evidently, in these cases, the defect enhances local field near it in the semiconductor, which increases the absorption dissipation.Even so, the LUSMP mode still suppress backscattering from defect.It should be indicated that in the lossless cases, the LUSMP can completely circumvent the defects II and III in the waveguide, as illustrated in Fig. 3d, f, h.We also numerically analyzed defects of different shapes, such as air columns with circular or triangle cross-sections, and the observed phenomenon is similar to that for air square columns.

Realistic waveguide using metamaterial
For the LUSMP waveguide, the UENZ material is a key component, whose relative permittivity is close to zero in the yz plane but larger than 10 in the x direction.Such a material does not exist in nature, however, it can be artificially realized using metamaterial technology.Here, as an example, the UENZ material is formed by alternating plasmonic and silicon layers, as depicted in the inset of Fig. 4b, and the period is much smaller than the wavelength.The plasmonic material may be a highly doped semiconductor, whose plasma frequency is significantly larger than that for the InSb.Thus, in the realistic waveguide, the UENZ layer comprises of seventy-five unit cells, and the period is 1 µ m.This waveguide can be fabricated with state-of-the-art technology: the thick Si film is first deposited on the substrate InSb, then thin plasmonic material and Si layers (of the metamaterial) are alternatively deposited, by using plasma-enhanced chemical vapor deposition (PEVCD) in a low-temperature environment 26 , and finally, a metal (e.g.Au) film is deposited on the top (thin) Si layer via conventional electronbeam evaporation 27 .According to the effective medium theory 28,29 , the UENZ metamaterial can be characterized by an effective (relative) permittivity tensor described by Eq. ( 1), and its elements have where δ is the filling ratio of the plasmonic material, ε p and ε q are the relative permittivities of the plasmonic material and silicon, respectively.To attain the desired values of ε z and ε x , we set δ = 0.09 , ε p = −116.8and ε q = 11.68.
We utilize the FEM to simulate wave transmission in this realistic waveguide.The UENZ thickness is taken to be 75 µ m (we chose a relatively small thickness to reduce the demand for computational resources), and the local model is considered.A magnetic current line source is placed in the middle of the UENZ layer to excite LUSMP in the system.The operation frequency is set at 2.12 THz.The simulated H-field amplitudes are plotted in Fig. 4a.We also numerically solved for modes in the realistic waveguide along the x-axis using FEM, and the results are shown in Fig. 4b, where the H-field profile is almost uniformly throughout the designed UENZ layer.A comparison with the corresponding homogeneous UENZ layer case, represented by a dashed line in Fig. 4b, (9)   shows good agreement.Therefore, the LUSMP in the realistic waveguide possesses a transverse modal size of approximately 75 µ m.In order to confirm the robustness of LUSMP in the realistic waveguide, a square silicon column with a sidelength of 15 µ m is introduced as a defect and placed at a distance of 150 µ m away from the source.The results are presented in Fig. 4c.It can be observed that LUSMP well circumvents the defect, and this is more clearly illustrated in Fig. 4d, which shows the distributions of H-field amplitude along the middle line of the UENZ layer in Fig. 4a,c.This confirms that LUSMP in the realistic waveguide is a robust unidirectional mode.

Applications of terahertz LUSMPs
We finally numerically demonstrate several advantages of LUSMPs over RUSMPs.Firstly, LUSMPs have a large modal area on the wavelength scale, enabling them to be efficiently excited by waves incident directly from air.
To show this, we use a Gaussian beam with a width of p to excite LUSMP in a waveguide with d 2 = 2 p , and the simulated H-field amplitudes are illustrated in Fig. 5a.It is found that 70.1% of the incident power is coupled to the LUSMP.In contrast, the RUSMP only has a coupling efficiency of 4% .Secondly, the modal size of LUSMP can be tailored by adjusting the UENZ thickness, as illustrated in Fig. 5b, where the modal size is abruptly varied from 2 to 0.2 p , and then back to 2 p .During the whole propagation process, the wave energy of LUSMP is first squeezed into a narrow channel and field amplitude is largely enhanced, then followed by almost full recovery with slight power loss due to the material loss ( ν = 0.002ω p ). Finally, LUSMP can be completely trapped with a hot spot of wavelength size.To achieve this, we have employed a backward-LUSMP waveguide to block the forward LUSMP, and simultaneously inserted a thin PMC slab with a thickness of 0.1 p between two InSb lay- ers under opposite magnetic fields to prevent wave leakage.The empty space between the UENZ layers in the two waveguides is filled with a dielectric with a relative permittivity of 18.The simulated H-field amplitudes are displayed in Fig. 5c, revealing the presence of a hot spot with a transverse size of 150 µ m.The large size of the hot spot avoids the risk of device burning-out.
Based on their large modal sizes and unidirectional propagation nature, LUSMPs can be utilized to develop a terahertz free-space isolator, enabling the transmission of plane waves in only one direction.The free-space isolator is comprised of two parts, as illustrated in Fig. 6.The lower part is an array of LUSMP waveguides of length 2 p , and the upper part is a dielectric layer of thickness 1.25 p .This dielectric layer has a very low relative permittivity of 0.004, and it is used to shape the output wavefront from the waveguide array.In each LUSMP waveguide, the InSb layer has a finite thickness of 0.05 p , which is far smaller than the UENZ-layer thickness of 2 p .For the proposed isolator, its frequency window is determined by the frequency range of the LUSMP waveguides.In the frequency window of this isolator, plane waves are prohibited to transmit through it when they are normally incident from its top, as illustrated by numerically simulated magnetic field distribution in Fig. 6a.On the contrary, plane waves are allowed to transmit through the structure when they are normally incident from below, as illustrated by numerically simulated magnetic field distribution in Fig. 6b. Figure 6c,d show the transmission efficiencies of the isolator for backward (incidence from top) and forward (incidence from below) propagation directions.At the central frequency 2.12 THz of the operation window [indicated by the shaded areas in Fig. 6c,d), the efficiency of forward transmission is nearly 88% .This also implies that incident plane waves can efficiently excite LUSMPs in the waveguide array, owing to their flat field profile.For backward transmission, our numerical calculations show that the isolation ratio reaches up to 60 dB at 2.12 THz.

Conclusion
In summary, we have conducted a comprehensive analysis of the guiding properties of a waveguide composed of metal, UENZ material, silicon, and magnetized InSb.It has been shown that in the upper bandgap of the magnetized InSb, the waveguide can support LUSMP mode, whose field is almost uniformly distributed over the thick UENZ layer.This LUSMP mode can possess a mode size larger than the wavelength.Moreover, the LUSMP mode is robust against the nonlocal effects of the materials, since the upper bandgap opened by the external magnetic field is topologically nontrivial.In comparison to (true) terahertz USMPs reported previously, which

Figure 1 .
Figure 1.(a) Schematic of the metal-UENZ-Si-InSb waveguide.(b) Dispersion diagram for the waveguide.Solid and dashed lines represent the dispersion curves for SMPs in the local ( β = 0 ) and nonlocal ( β = 1.07 × 10 6 m/s) models, respectively.The dotted lines are light lines in the UENZ material, and the green shaded areas represent the bulk-mode zones in the magnetized InSb for the nonlocal model.But the lower limit of the upper bulk-mode bandgap for the local model is indicated by the dot-dashed line.In the nonlocal model, the lower bulk-mode zone covers the upper bandgap, thus it becomes an incomplete bandgap.Moreover, USMP vanishes in the lower bandgap.The parameters of the Si and UENZ layers are ε z = 0.001 , ε x = 50 , d 1 = 0.63 µ m, and d 2 = 300 µ m, and the plasma frequency of the InSb is f p = ω p /2π = 2 THz.(c) Magnified view of the USMP window in (b).The unidirectional window in (c) ranges from 2.06 THz to 2.2 THz.(d) Dispersion curves for SMPs in the local (solid) and nonlocal (dashed) models, where ε x = 300 and d 1 = 0.18 µ m, and the other parameters are the same as in (b) and (c).In this case, the whole upper bandgap is a unidirectional window.

Figure 2 .
Figure 2. (a) H-field profile of LUSMP for d 2 = 300 µ m.As a leaky mode, the H-field amplitude of USMP in the nonlocal and lossless case does not decay to zero in the InSb layer.(b, c) Dependence of k i on d 2 in the nonlocal ( β = 1.07 × 10 6 m/s) and local ( β = 0 ) models for various ν values.The frequency is all 2.12 THz in (a), (b) and (c).(d) Propagation lengths ( L = 1/2k i ) of LUSMP and RUSMP as functions of frequency.Various ν values are considered.For the LUSMP waveguide, the UENZ layer has d 2 = 300 µ m.The other parameters are ε z = 0.001 , ε x = 50 , and d 1 = 0.63 µ m for all the cases.

Figure 3 .
Figure 3. Simulations for waveguides in the nonlocal model.(a-f) Simulated magnetic field amplitudes launched by a magnetic current line source at f = 2.12 THz for different cases.(a) Lossy case without defect, (b) lossy case with defect I, which is a square air column with side length 75 µ m, centered in the middle of the UENZ layer.Note that the field patterns for the lossless cases without defect or with defect I are almost identical to those in (a) and (b).(c, d) Lossy and lossless cases with defect II, which is the defect I shifted down by 75 µ m. (e, f) Lossy and lossless cases with defect III, which is the doubled defect I. (g, h) Distributions of magnetic field amplitude along the middle line of the UENZ layer for the lossless and lossy cases without or with defects.For all the lossy cases, ν = 0.002ω p .The other parameters are the same as in Fig. 1b.

Figure 4 .
Figure 4. (a, c) Simulated H-field amplitudes in the realistic waveguide without and with defect.The UENZ layer in the waveguide consists of alternative plasmonic and Si layers, and the period is 1 µ m.The defect is a square Si column of side length 15 µ m.(b) H-field distribution along the x axis.The dashed line represents the corresponding H-field distribution for the effective model where the UENZ layer is homogenized.The inset is a schematic of the realistic waveguide structure.(d) Distributions of H-field amplitudes along the waveguide axis in (a) and (c).The frequency is 2.12 THz for all cases.

Figure 5 .
Figure 5. Simulated H-field amplitudes in systems where LUSMPs are coupled or manipulated at f = 2.12 THz.(a) LUSMP is efficiently coupled by an incident Gaussian beam with width p .(b) Modal size of LUSMP is tailored by varying the UENZ-layer thickness.It is sharply reduced from 2 p to 0.2 p , and then back to 2 p .In the middle section, the modal energy is largely squeezed.(c) Complete trapping of LUSMP.The transverse size of the generated hot spot is equal to p .The loss of the semiconductor is ν = 0.002ω p , and the other parameters are the same as in Fig. 1b.

Figure 6 .
Figure 6.Terahertz free-space isolator formed by an array of LUSMP waveguides.(a, b) Simulated magnetic field amplitudes for plane waves incident from the top and bottom of the isolator.(c, d) Transmission efficiencies of the isolator for the forward and backward directions.The loss of the semiconductor is ν = 0.002ω p , and the parameters are the same as in Fig. 1b.